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The Lund string fragmentation model is applied in a non-perturbative calculation of the invariant 
production cross sections of pions from proton-proton collisions in the soft pr region. Invariant 
production cross sections of pions and kaons from proton-proton collisions in the hard pr region are 
calculated from the Feynman-Field perturbative QCD parton model. Parameterizations of these 
invariant production cross sections are described. 
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q , I. INTRODUCTION 

. This work was originally motivated by the need of parameterized meson production cross sections of pp — ► hX 
reactions (p for primary proton, h for identified hadron production and X for unidentified hadron production) for a 
NASA nuclear transport code called HZETRN The parameterization scheme presented here parameterizes theory. 
Experimental data are used mostly as checks. The parameterizations are based on two main calculations — a non- 
, perturbative QCD string fragmentation Lund model calculation in the soft px region (px < 1 GeV) and a perturbative 
QCD Feynman-Field parton model calculation in the hard pr region (pr > 1 GeV). The threshold of 1 GeV separating 
the soft and hard pt regions is chosen to be the proton mass m p — 938 MeV. Descriptions of the Lund and Feynman- 
Field models can be found in original sources |7j and a review paper [0 . Both the Lund model and the Feynman-Field 
model are stochastic models and are not explicitly quantized. In this sense, both are phenomcnological models. The 
goal of this work is to calculate E d 3 a/dp 3 with sufficient accuracy using reasonable theoretical models so that reliable 
— ■«. " parameterizations of cross section formulas can be determined. The Feynman-Field model is implemented numerically 
'-^ \ by the Monte Carlo integration package VEGAS |fl). The invariant cross sections are parameterized for pions in 
II ■ the soft pt region and for pions and kaons in the hard pt regions. 
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String fragmentation models such as the Lund Model fit experimental data well. String theory reproduces the 
linear potential predicted by non-perturbative QCD as in lattice gauge field theory. These observations hint at the 
possibility that QCD string may be conducive to solving non-perturbative QCD. Although the Lund model is a (1 + 1) 
model, it reproduces the essential dynamics of the system as long as information on angular momentum is not needed. 
Typically the Lund model is implemented numerically using Monte Carlo simulation in JETSET and PYTHIA ||. This 
section shows how to calculate meson production cross section formulas analytically in the non-perturbative QCD 
region using the Lund model. Special attention is given to the invariant production cross sections of pions from 
proton-proton scattering in the soft pt region (px < 1 GeV) where non-perturbative effects dominate. The results of 
the string model are compared against inclusive pion production cross section data of proton-proton collision. 

This section uses the same notations of Reference [[| . The basic result of the Lund Model is the "area law" which 
is summarized as (|] 

j n / n \ 

dP = ds _ (i _ z y Y[Ndp 2 0j 5+(p 2 0j -m 2 )sl Prest - Y^Poi e~ bA . (1) 
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P is the probability of the string fragmentation process, b and N are constants and s is the total energy square of n 
produced mesons. Lightcone coordinates are used in the Lund model, x± = t ± x and p± = E ±p. The momentum 
Prest = (W-, W+) is that of the parent virtual quark pair and poj is that of the j-th rank meson. Momentum transfer 



is 



-t = —q = W-W+ where t is a Mandelstam variable. The mass of the system is m. The area A = T + A r 



of the polygon in Fig. [I] has a geometric interpretation as the residual energy of the virtual quark pair following the 
fragmentation process. T — k 2 x+x~, A and A res t are Lorentz invariant kinematic variables. The symbol z — ■ z j 
denotes the sum of the fractions zoj of lightcone energies of the produced mesons along the x+ direction. The quantity 
r defines the surface of constant proper time along which the string is broken. The area law captured in Eq. (|l|) is 
usually incorporated into Monte Carlo simulation programs such as JETSET and PYTHIA ||. Monte Carlo subroutines 
are too slow for nuclear transport codes. Computational constraints motivate the search for parameterized production 
cross section formulas. 

By taking dP = da/ao, where <7o is the total cross section, and using the identities found in References ||, ||, 
Eq. ([!]) can be rewiritten as 



da = er ds dz (1 - z) a e br 5 \z - z oj \ 8 \ s — 



-oj 



n , 

| J A — — ^ e- bA " 



(2) 
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where z j is the fraction of the parent quark momentum p+ that goes into the momentum p j+ of j-th rank meson. 
For n > 1, the area A rest can be divided into n rectangles so that 



Arest — ^ ^-oj • 



(3) 



3 = 1 



Fig. illustrates the first of such n rectangles and the geometric properties of the measures of its sides. By utilizing 
the relation 



Poj + 



Poj- 



it can be shown that 



A Q j z j r. 



When n = 1, there are 2 vertices and A rest — z a iT. In this case, Eq. (H) can be integrated over s and z to give 



da = (To N (1 - z) a e 



a -b(l+z)T 



(4) 



(5) 



(6) 



It is easy to show that the Feynman variable x = p z /{p z )max is equal to z when there is only 1 produced meson so 
that dx = dz. The differential cross section is found to be 



da 

= (Jo N < 

dx z 



( 1 ~ Z T p -b(l+z)l 



The invariant quantity 



r = m h e^y (W± - m h e ±y ) , 



(7) 



(8) 



where rrih is the mass the produced meson, can be understood intuitively through geometrical considerations by 
focusing on the rectangle in the spacetime diagram next to the first rank vertex V\ in Fig. [|. Generally, the lightcone 
energies of the parent quarks can be taken to be the momentum transfer W- — W+ — q. The feynman variable can 
be defined as 



p z 2to/,. sinh y 

(Pz)max q 

With Eqs. (||) and (g), Eq. (0) can be rewritten as 

(j; ~ p -b{l+z)m h e-«(q-m h e"") 



da 

dx 



a N' 



z 



,-b(l+z)q(q-mh e y )a 



ao N Q ^ e -Hl+*)mh(q-m h ) 

z 



,-b(l+z)q(q-mh)x 



Ce 



-Dx 



(9) 
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Eq. ( pH ) is obtained by setting y — 0. In general, experimental data average over rapidity that typically centers 
around zero such that y — is a reasonable simplification. In addition, average values of internal variables z and q 
are used in Eq. ( |TT| ) so that C and D can be treated as constants when comparing with experimental data. The form 
of Eq. (|lj) is consistent with experimental data as shown in Figs. § and |[ 

When n = 2, Eq. (0) is transformed as follow by integrating over z Q \ and z Q i and letting £ = z Q i(z — z Q i): 



der = (To iV (1 — z) a dz 5(z — z \ — z 2) ds S [ s — 

x e ~br e ~b( Z01 +z o2 )r 

Zol z o2 

= (j a N (I- z) a dzdsd 
= a N(l- z) a dzdsS 
= (T N (1 - z) a dz ds - 



Zol 



Zo2 



Zol 



Zol 
d£ 



dzoi 



Z ol (z - Z l) 
_ p -b(l + z)T 



-b(l+z)T 



^ J ^z 



,-6(i+z)r 



7 2 _ 4m 2 



(13) 



One of the typical assumptions of the form of momentum transfer in QCD is —t 
s ~ -t - (1 + z)T = Ap^, - (1 + z)F, we obtain 

ds = idp^, 

by keeping z and T constant. On the lightcone, p z = E so that 

dE dp z 



dz = 



E 



E 



Eqs. Q and (|l|) together give the relation 



dz ds 



dp\ dp z 



dp 3 



E E 

By combining Eqs. (|l^) and (|l6|), the invariant production cross section can finally be shown to be 

E fa = 4iV(l-z) a (To e ~ H i +z) r^ 
dp 3 y / s 2 z 2 - Asm 2 

If W+ in Eq. (^) is taken to be q — 2px, Eq. (JT^) can be expressed in terms of pr as 

d 3 a 



E- 



dp 3 



AN{l-z) a a Q cHl+z)m i 



v s^z 



2 r 2 



As m 2 



= At 



4N(l-z)°<r f>b(1+a)< 
\J s 2 z 2 — As m 2 

-Bp T 



5 -2b(l+z) m h e-"p T 



-2b(l+z) rrih Pt 



Ap\. With the relations 



(14) 
(15) 

(16) 

(17) 

(18) 

(19) 
(20) 



Again y = is assumed in Eq. (|T^) . Average values of s and z are used when comparing with experimental data so 
that A and B are treated as constants. It must be emphasized that s is not the beam energy square but the total 
energy square of the produced mesons in the string fragmentation process. The parameters for pion cross section are 
listed in Table |. B is extracted from data as shown in Figs. |]-|6|. A is chosen to match the pion cross sections at the 
boundary of pr — 1 GeV between the soft and hard pr regions so that 



A = 3x l(T 28 e s 



(21) 



in unit of cm 2 /GeV 2 for pions. These figures suggest that the cross sections are approximately constant in yfs and 
9 cm in the range of energy being parameterized (22 GeV < s/s < 63 GeV). Therefore the present parameterization 
of pion production cross sections for px < 1 GeV is expected to be valid for the energy range not narrower than 
20 GeV < y/s < 70 GeV. Kaons are not analyzed nor parameterized because of insuffucient data. 
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III. PION AND KAON PRODUCTION IN THE HARD p T REGION 



The invariant cross sections of pions and kaons from proton-proton scattering in the hard px region (j?t > 1 GeV) 
are treated in this section by using perturbative QCD. The theoretical foundation of the present calculation is based 
on the Feynman-Field parton model. A comprehensive review of the model is given in References @, §. The cross 
sections are parameterized at the end. 

The Feynman-Field invariant production cross section formula incorporates the parton distribution functions, 
fA/a(x a ,Q 2 ), obtained from DIS experiments and the fragmentation functions, D^(z,Q 2 ), derived from a combi- 
nation of stochastic arguments and parameterizations of data. The cross section formula is given as M 



ap a it *■ — ' 



a , b 



dx a 



dx b f A /a(x a , Q 2 ) fB/b{xb, Q 2 ) D^(z c , Q 2 ) — 

z c dt 



(22) 



with 



xi 



1 — X2 ' 

XgX 2 



(23) 
(24) 



Monte Carlo integration package VEGAS is used to calculate Eq. ( [22] ) . The parton distributions of proton is given by 
the CTEQ6 package |l3). The QCD running coupling constant, cx a (Qr), is the renormalized coupling constant described 
in Reference ||. A typical value of A = 0.4 GeV is used inside a s {Q 2 ). The internal scattering cross sections of the 
reactions — ► gg and gg — > gg are excluded from the integral because gluons do not fragment into hadrons. The 
fragmentation functions used for this calculation are the original fragmentation functions of Feynman and Field [jl4| . 
For the pp — > n X reactions, the fragmentation functions are 



DZ (z) = Dl (z) 



(n + l)(l-z)", 



(25) 



Df(z)=p 2 



(n + l)(l-zY 



(26) 



Df{z)=DZ (z) 



P + (3 2 



1 - z 



(n+1) (!-«)", 



Df(z)=Df (z) = DZ~ (z) = 



2 1 1 ~ ' 1 (n + l)(l-z)", 



and for pp — > K X reactions, the fragmentation functions are 



D« + {z)=Df{z) 



1 



P + P 2 



1 - z 



(n + 1)(1- z) n , 



(27) 
(28) 

(29) 



D? \z) = Df(z) = Df{z) = Df(z) = L ' ' ; 



£»f {z) = D*-{z) = 



(n+l)(l-zr, 



(n+l)(l-zY 



DT (z) = DfT (z) = Df (z) = Df (z) = £ H—^j (n + 1) (1 - zf 

where (3 = 0.4. The distributions of c, b and t quarks are sufficiently low that 

D h c (z) = DZ(z) = D?(z)=0, 



(30) 

(31) 
(32) 

(33) 
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for any hadron h. Feynman fixed n — 2 in his original paper. In this work, n is a parameter freely adjusted to fit data. 
There is a subtlety involved in summing all the parton contributions over a and b in Eq. (|2^) that is related to the 
relative probabilitistic nature of the parton distributions and our ignorance of the number of sea quarks and gluons 
inside the proton. The parton distributions are normalized to unity so that they give only the relative distributions of 
the partons. The parton distributions give only the ratios of the partons in a hadron but not their numbers. In order 
to sum over a and b partons in Eq. (g2J) correctly, an integral multiplicative constant for each of the hadrons A and 
B must be provided. These multiplicative integral constants are not known a priori but are determined a posteriori 
by fitting data. In other words, Eq. (|22| ) can be modified as 



E- 



dp 3 



N A N B 



E 

{a,b} 



dx„ 



dxh 



x fA/a( x a,Q 2 ) f B /b(Xb, Q 2 ) D h c {z Cl Q 2 ) 



1 da 

z c dt' 



(34) 



where Na and Nb are the multiplicative constants corresponding to fA/a(x a ,Q 2 ) and fB/b{xb,Q 2 ) respectively and 
S{ a b} ^ s t ne sum over the parton types instead of a sum over the partons per se. If A = B, the overall multiplicative 
constant, Na Nb, is an integer square. If A ^ B, the overall constant, Na Nb, is still an integer. In the case of fitting 
pQCD calculations to experimental data of a pp —> ttX reaction, a factor of 100 is missing if one simply sums over the 
parton types. It implies that the multiplicative constant for the parton distributions of proton is N p — 10. The cross 
sections for K production is approximately half of that of n production. It implies that the multiplicative constant 
for the pp — > K X reaction may be N p = 7. For the purpose parameterizing the shape of the kaon production cross 
section, an exact scale is not required. Therefore N p in the pp — > KX reactions is arbitrarily set to be the same as 
that of pp — ► ttX reactions at N p = 10. This choice is adequate because fits to experimental data of kaons are not 
being pursued in this work due to the lack of experimental data for kaons. 

It is observed that the invariant production cross sections have the same basic shape regardless of the reactions, 
i.e. an exponentially decaying function of the form exp(-ax^) at low pt and a suppression at high px which drops 
off to zero before the edge of suppresion at px < y/s/2. The cross section is at its maximum at pr = and decreases 
monotonically in pt- The Feynman-Field code used in this calculation assumes that Q 2 = <^p\. In other words, by 
combining the previous two statements, the cross section is at its maximum at Q 2 — and decreases monotonically in 
Q 2 . This observation indicates that hadron fragmentation is more favorable at low Q 2 in that a parton preserves more 
kinetic energy to be made available for hadron fragmentation. It is also observed that the cross section is suppressed 
at high pt and that the edge of suppression of the cross section is at low pt and gradually increases toward 

high pt- The reason for this phenomenon is mostly due to the choice of Q 2 = 4px such that s > Q 2 or equivalcntly 
%/s/2 > pt- It turns out that the edge of suppression along pt is a function of yfs. The comments made so far apply 
to any angle 6 cm when Js is replaced by sin The basic features of the graphs of the invariant cross section is 
discussed in Reference g . 

A lot of effort has been invested in parameterizing E d 3 a/dp 3 from experimental data for use in the HZETRN code [[l2| . 
This work takes a different approach by parameterizing theory. Experimental data are used merely as a means to fine- 
tune the parameterizations. Monte Carlo integration is the fastest numerical integration scheme available but it is not 
faster than an analytic formula. On the other hand, the explicit computation of the double integrals in the Feynman- 
Field model is a daunting task if tractable at all. These constraints motivate the present parameterization scheme. 
The method of finding the parameters is mostly one of trial-and-error. Guesses of the appropriate functions are made 
for different parts of the curve. The pieces are put back together at the end and refit until the parameterization has 
the desired global properties. The final form of paramterized cross section formula is found to be 



PT, 



= ^4 e -35.4/3(p£-l) 







( 




X 


exp 


1 — CSC 













Pt ~ 1 



(sin9 crn ^/s -i 




(35) 



for pt — 1 < a 



sin 9 C m \/s 



— 1 ] . It is taken that 



dp 3 



(36) 
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for pt — 1 > a ( sm 6c ™ ^ — ij . A is a scale factor such that 



A = E ^-g-Os/i, p T = 1, #cm - 90°). (37) 



dp 3 

The functions a and /3 are 

a = 2-e~ Qsin9 -^, (38) 



/9 = #>+A (sin0 cro Vs) *, (39) 

and the parameters a, b, /3o, (3\ and /32 are freely adjusted to fit the curve. The unit of energy- momentum is GeV 
for the present parameterization. A is determined by varying n to fit the data with the Feynman-Field calculation. 
The first exponential in Eq. ( |35|) controls the shape of the curve at low px and the square bracket term controls the 
suppression at high px- The function a shifts the edge of the suppression so that the edge is located at yfs/2 at low px 
and y/s at high px- There is a threshold px ~ 0.2 GeV set by the CTEQ6 package. In addition, non-perturbative effects 
become more prominent in the soft px region (px < 1 GeV) so that the pQCD code cannot be applied there. For 
these reasons, the present parameterization focuses on the region px > 1 GeV. The rapidity distributions of hadron 
production is typically symmetric around y = 0. From Eq] 



r} = - ln[tan(0/2)] = - In 



m T cosh y — m? + mx sinh y 



m T cosh 2 y — m 2 — mx sinh y 



(40) 



Vf^l. (41) 



and 

V= 2 ~ \E- Pz 

it can be easily shown that y = 0, Xl = Pz/iPz) max — and 9 cm — 90° are the same statements. Many experiments 
average the data over a range of y or xl symmetric around zero. In these cases, the average center-of- momentum angle, 
9 cm , would be 90°, which happens to be the most prominent contribution according to Fig. 0. Hadron productions 
in space radiation problems are mostly in the forward direction or equivalently 9 = 0. In highly relativistic regimes, 
9 cm — 90° is transformed to a small 9i a b- For example, 9 cm = 90° is equivalent to 9i ao = 9.17° at yfs = 70 GeV [fT6| . 
The scattering angle 9i a \, decreases as yfs increases. At sufficiently high y/s, 9i a ^ is effectively zero so that a one- 
dimensional transport code is justified and that the largest angular contribution of the cross section is that at 9 cm = 
90°. 

Figs. ||-jl3|show the goodness of fit between the parameterized cross sections and the pQCD calculations at 9 crn = 90° 
and Fig. |Ij shows the goodness of fit at various angles. Apparently the same parameterization of j3 works for all 
types of pions and kaons (see Figs. ||-|l3|). The Feynman-Field model fits data well when n is allowed to vary as it 
is illustrated in Fig. [l5|. Generally speaking, the parameterization of the model at a fixed value of n = 2 does not 
always fit data well. For this reason, data are parameterized separately by reparameterizing (3. The parameterizations 
of pions are illustrated in Figs. p~6|— |l8|. Data of kaons are fragmented so that only the theory (but not the data) is 
parameterized. The parameters are tabulated in Table [ll| The parameterization of data is distinguished from that 
of theory by labelling the paramters of the former as f3 and those of the latter 0' . All other parameters are the same 
for both data and theory. A sample graph of the fit between the present parameterization and data with angular 
dependence is shown in Fig [19[ 

The parameters in Table |II| are obtained by fitting curves by hand. Curve-fitting algorithm such as the Levenberg- 
Marquardt does not work well because of the presence of singular matrices inherent in the present model. At this point, 
there is not a cleverer method to find these parameters automatically. When the parameters are obtained by hand, 
several data sets of E d 3 a/dp 3 versus px at different energies with px > 1 GeV and 9 cm = 90° are needed. In general, 
more data sets means smaller \ 2 - I n the case of parameterizing theory, the Monte Carlo program can generate as 
many theoretical data sets as needed. In the case of parameterizing experiment, data are generally fragmented except 
those of pions. At least 3 experimental data sets are needed to fit (3. Pion data are generally quite copious. In the case 
of kaons, the scarcity of data prevents the parameterization of their experimental fits. The present parameterization 
of pion production cross sections for px > 1 GeV is expected to be valid over 10 GeV < sin9 cm y/s < 70 GeV, which 
is given by the range of energy of the experimental data used in this analysis. 

Most experiments agree with the shape of theoretical curves. However there are diagreements in the magnitudes 
of the cross sections among experimental data sets, sometimes even by the same authors. Figs. O and show that 
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Abramov et al. published 2 sets of pion data in identical energy regimes in 2 consecutive years that are different by 
3 orders of magnitude [jl6|, |i~7| . This phenomenon occurs quite regularly in Ed 3 a/dp 3 data. Some discretion must be 
exercised in choosing the experimental data to parameterize. 



IV. CONCLUSION 



This work shows how to use the Lund model to calculate the invariant production cross section non-perturbatively 
in the soft px region. This model predicts that the functional form of the cross section in this sector is a simple 
exponential. This prediction is confirmed by experiment in the case of pions for 8 = 90°. The cross section formula in 
Eq. (|2^) has no angular nor energy dependence. Although the prediction of angular independence is not yet confirmed 
by data, energy independence is shown to be approximately correct according to the Figs. ^j^. The Feynman-Field 
model is a work horse for calculating inclusive cross sections and is generally accepted to be an accurate description 
of data in the high px region. Several typical assumptions of pQCD have been adopted in the present calculation 
that are not necessarily unique, such as the forms of Q 2 and the fragmentation functions. The form of momentum 
transfer is assumed to be Q 2 = Ap T in the present work. In principle, other guesses, such as 

^ 9 2stu , , . 

Q = N-7i > ( 42 ) 
s z + t l + u z 

can also used. Nevertheless it is unlikely that the use of Eq. ([42]) in the Feynman-Field model will drastically change its 
predictions. More experimental data of kaons are needed to determine the corresponding paramterizations. Otherwise, 
the parameterization procedure can in principle be applied to all types of mesons. The production cross sections of 
baryons arc also needed for a realistic nuclear code. It is not yet clear that the Lund model and Feynman-Field can 
be easily modified to incorporate baryon production. Eventually production cross sections of hadrons from pA and 
AA collisions must also be included in the transport code. The phenomena of heavy nuclei collisions are much more 
complicated because many-body effects such as the EMC effect, Cronin effect, nuclear shadowing, jet quenching and 
gluon-plasma phase transition have to be considered. It is hoped that the present work provide some basic ideas for 
more complicated calculations to be undertaken in the future. 
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FIG. 1: A spacetime diagram of the (1 + 1) Lund string fragmentation process. The vertices Vi denote the spacetime points of 
the creation of virtual quark pairs. The quarks are massless so that their trajectories are light-like. The arrows indicate the 
directions of the trajectories of the produced hadrons. 
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FIG. 2: Comparison between experimental da/dx and the string model result. The constants used in the exponential function 
are chosen to fit the data and are not parameterizations. The data are published in Reference M. 




FIG. 3: Comparison between experimental Eda/dpL and the string model result. The constants used in the exponential 
function are chosen to fit the data and are not parameterizations. The data are published in Reference J4[. 




FIG. 4: Comparison between experiment and theory of Ed 3 a /dp 3 for pp —* tv + X in the soft pr region (pr < lGeV). The data 
are published in Reference pi. 



• V~s = 22 GeV 
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FIG. 5: Comparison between experiment and theory of E d 3 a/dp 3 for pp —> ty~ X in the soft pr region (pr < lGeV). The data 
are published in Reference pi. 
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FIG. 6: Comparison between experiment and theory of E d 3 a/dp 3 for pp — > ty° X in the soft pr region (pr < 1 GeV). The data 
are published in Reference pi. 
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FIG. 8: Parameterization of the Feynman-Field model for pp — > ir + X at 6 cm = 90° and n = 2. The values of listed in the 
figure match the curves from left to right. 
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FIG. 9: Parameterization of the Feynman-Field model for pp — > 7r~ X at 6 cm — 90° and n — 2. The values of listed in the 
figure match the curves from left to right. 
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FIG. 10: Parameterization of the Feynman-Field model for pp — > tv°X at 9 cm = 90° and n — 2. The values of ^fs listed in the 
figure match the curves from left to right. 
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FIG. 11: Parameterization of the Feynman-Field model for pp — > K + X at 9 cm = 90° and n = 2. The values of listed 
the figure match the curves from left to right. 
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FIG. 12: Parameterization of the Feynman-Field model for pp — > K°X at 6 cm — 90° and n = 2. The values of listed in the 
figure match the curves from left to right. 
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FIG. 13: Parameterization of the Feynman-Field model for the reactions pp — > K~ X and pp — > 'x at 6 cm = 90° and n = 2. 
The values of -/s listed in the figure match the curves from left to right. K~ and K have the same fragmentation functions 
and hence the same Feynman-Field cross section. However their experimental cross sections may not be the same. 
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FIG. 14: Parameterization of the Feynman-Field model for pp — > tv + X at various angles and n = 2. The values of the angles 
listed in the figures match the curves from left to right. 



le-27 



le-30 - 



% 



Oh 

-a 



^ le-33 
pq 



le-36 - 



n V. 



pp— >n X 



\ v. 



e 



"OS* 



:90° 



□ 
A 



= 1 1.5 GeV (Abramov 80) 
V~s = 11.5 GeV (Abramov 81) 

= 38.7 GeV (Jaffe 89) 
VT= 52.7 GeV (Busser 76) 
FF model, variable n to fit data 
FF model, n = 2, VT= 11.5 GeV 
FF model, n = 2, = 52.7 GeV 
FF model, n = 3, -fs = 52.7 GeV 



> n = 0.5 
V 

\ 

\ 



J !, L 



n = 4 



J i L 



. N L 



5 

P T (GeV) 



10 



FIG. 15: Comparisons of the Feynman-Field model fit for pp - 
sets Abromov80, Abramov81, Jaffe89 and Busser 76 are Jig], [117 



7r + X for various n at 6 cm 
, |Q and H respectively. 
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FIG. 16: Comparisons of pQCD and parameterized cross sections with pp — * ir X experimental data at 6 C m = 90°. The 
parameter n in fragmentation functions of the Feynman-Field model are adjusted freely as shown in the graph to fit the data. 
The references of the data sets Adams96, Demarzo87 and Akesson89 are Eq] , J2l| and [pi) respectively. 
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FIG. 17: Comparisons of pQCD and parameterized cross sections with pp — > ix + X experimental data at 9 cm — 90°. The 
parameter n in fragmentation functions of the Feynman-Field model are adjusted freely as shown in the graph to fit the data. 
The references of the data sets Abromov80, Abramov81, Jaffe89 and Busser76 are jlij , |l7| , flit and [^| respectively. 
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FIG. 18: Comparisons of pQCD and parameterized cross sections with pp — > n~ X experimental data at 8 cm = 90°. The 
parameter n in fragmentation functions of the Feynman-Field model are adjusted freely as shown in the graph to fit the data. 
The references of the data sets Abromov80, Abramov81, Jaffe89 and Busser76 are jlij , |l7| , flit and [^| respectively. 
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FIG. 19: Comparisons pQCD and parameterized cross sections with pp — > ty + X experimental data at various energies and 
angles. The parameter n in fragmentation functions of the Feynman-Ficld model are adjusted freely as shown in the graph to 
fit the data. The experimental data set is published in Reference [E3[. 
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TABLE I: Parameters of the parameterized cross section formula in the soft pr region [pr < 1 GeV). A is chosen to match the 
pion cross sections at the boundary between the soft and hard pr regions and has the unit cm 2 /GeV 2 . B is extracted from 
data as shown in Figs 





7T° 


7T + 


7T 


A 


4.45E-26 


6.64E-26 


6.64E-26 


B 


5.0 


5.4 


5.4 
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TABLE II: Parameters of the parameterized cross section formula. A is chosen to fit data and has the unit cm /GeV 2 . ft 
parameterizes data while ft parameterizes the Feynman-Field model at n = 2. ft of kaons are not parameterized because of 
insufficient data. 







TV 


+ 


7T~ 


K 


K 


K~ 


IF 
K 


nf 

A 


0.33 


0.33 


0.33 


0.33 


0.33 


0.33 


0.33 


at 

ft 


1.9339 


1.9339 


1.9339 


1.9339 


1.9339 


1.7931 


1.7931 


02 


1.0558 


1.0558 


1.0558 


1.0558 


1.0558 


0.9849 


0.9849 


a 


4.855E-3 


4.855E-3 


4.855E-3 


4.855E-3 


4.855E-3 


7.5E-3 


7.5E-3 


b 


0.98 


1.00 


0.90 


1.00 


0.90 


0.85 


0.85 


A 


3c-28 


3e-28 


3e-28 










ft 


0.30 


0.20 


0.20 










Ih 


0.3337 


0.3228 


0.3510 










A 


0.3774 


0.1472 


0.1815 











